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Abstract.  Barrier  calculations  based  upon  solutions  of  the  Liouville  equation  in 
the  coordinate  representation  reveal  a  complicated  spatial  dependence  of  the 
quantum  distribution  function  near  and  within  the  barriers.  Within  the  framework 
of  classical  transport  this  spatial  dependence  suggests  equilibrium  electron 
temperature  values  that  differ  from  the  ambient.  The  prospect  of  quantum  heating 
and  cooling  under  equilibrium  conditions  is  examined  and  disF>elled  in  favour  of 
an  interpretation  that  includes  density-gradient  contributions. 


1.  Introduction 

Calculations  based  upon  solutions  of  the  Liouville 
equation  in  a  density  matrix  formulation  yield  a  com¬ 
plicated  spatial  distribution  with  a  mean  kinetic  energy 
in  equilibrium  that  differs  significantly  from  the  classical 
result.  In  particular,  where  classical  physics  teaches  that 
the  energy  per  particle  is  kg  r/2  per  degree  of  freedom  for 
a  Boltzmann  distribution,  quantum  physics,  as  pointed 
out  by  Wigner  [1],  teaches  otherwise.  The  origin  of  this 
difference  lies  in  the  presence  of  quantum  mechanical 
forces  arising  from  gradients  in  density  (Ancona  and 
lafrate,  [2]),  and  are  suggestive  of  a  spatially  dependent 
local  temperature  in  both  equilibrium  and  non- 
equilibrium  cases,  although  spatial  dependent  carrier 
temperature  in  equilibrium  introduces  interpretive  dif¬ 
ficulties.  To  avoid  this  difficulty  one  either  abandons  the 
spatial-dependent  temperature  description,  or  retains  it 
for  non-equilibrium  studies  and  seeks  another  descrip¬ 
tion  for  equilibrium.  But  in  either  case,  it  is  necessary  to 
demonstrate  its  origin.  This  is  provided  below  for 
equilibrium  conditions. 


2.  Energy  and  temperature 

Qassical  physics  indicates  that  the  mean  kinetic  energy 
in  equilibrium  for  carriers  obeying  Boltzmann  statistics  is 

<£>=  J  d^p(p^/2m)f(x,  p)  =  Ipfeg  T  (1) 

where  /(x,  p)  is  the  classical  distribution  function,  and  e, 
the  mean  kinetic  energy  per  particle,  is  independent  of 
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position,  as  is  the  consequent  electron  temperature.  For 
quantum  structures,  in  which  quantum  distribution 
functions  are  required,  e  is  generally  spatially  dependent 
[1],  and  on  the  basis  of  equation  (1)  suggests  a  spatially 
dependent  equilibrium  carrier  temperature.  Because  of 
the  significance  of  carrier  temperature  in  interpreting 
hot-carrier  phenomena,  the  spatial  dependence  of  the 
mean  energy  per  particle,  and  the  origin  of  this  de¬ 
pendence  is  discussed  through  solutions  to  the  Liouville 
equation 

ihdpjdt  =  [//,  (2) 

which  in  the  coordinate  representation  is  a  differential 
equation  for  p(x,  x\  t) 

dp/dt  +  (fi/2miXV^  -  V^)p 

-(l/ifi)[K(x,t)- K(x'.t)]p  =  0.  (3) 

To  expose  the  essential  features  of  this  discussion,  we 
assume  Boltzmann  statistics,  spatial  variations  only 
along  the  x  direction,  and  free  particle  behaviour  along 
the  y  and  z  directions.  Transforming  equation  (3)  to 
centre  of  mass  and  non-local  coordinates,  r={x  +  x')/2, 
C=(x-x')/2,  p=>p{r  +  Q,r~Q,  we  find 

p.+W2mi)p^-{\/ih)lV{r  +  i;,  t)  -  K(r-C,  t)]p  =  0.  (4) 

In  the  above  equation  subscripts  denote  differentiation. 
The  potential  V  in  equations  (3)  and  (4)  is  the  sum  of  all 
heterostructurc  contributions,  Fo(x),  and  contributions 
from  Poisson’s  equation: 

d/dxl(^x)dVJdx]  =  -  <?^[p(x,  x)  -  Po(x)]  (5) 

where  the  subscript  ‘sc’  denotes  self-consistent;  Po(x)  is 
the  background  ‘jellium’  doping  distribution.  The  dia¬ 
gonal  components  of  solutions  to  equation  (4)  (along  the 
diagonal  r  =  x  and  C  =  0)  provide  the  density,  while  the 
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expectation  value  of  energy  <£>  is  obtained  from  the 
diagonal  components  of  the  kinetic  energy  density  matrix 
[3] 

£{r  +  r  -  ()  =  -  h^lSmp.^.  (6) 

An  approximate  form  of  the  expectation  values  of  the 
density  and  the  energy  density  [1,4,5]  for  one  degree  of 
freedom  is: 

p(x)  =  p(x,  X)  =  po  exp[ - (K  +  Q/3)/kB  7]  (?) 

E(x)  =  [(fee  T/2) + (A  V24mkB  T)  K„]p(x).  (8) 

In  equation  (7),  po  is  a  reference  density,  and  Q(x)  is  the 
Bohm  quantum  potential  (see,  e.g.,  [6]): 

Q(x)^ -(f>^/2m)(p'^\Jp''\  (9) 

The  s«:ond  term  of  equation  (8)  is  referred  to  as  the 
Wigner  contribution.  In  equilibrium  the  spatial  de¬ 
pendence  of  the  energy  per  particle,  c,  as  given  by 
equation  (8)  is  second  order  in  A.  To  this  order,  if  the 
potential  appearing  in  equation  (8)  is  represented  by  the 
Boltzmann  relation  between  density  and  potential 
energy,  p(x)=Poexp[—y(x)/kBT],  it  is  seen  that  the 
spatial  dependence  of  £  is  a  direct  consequence  of  the 
spatial  derivatives  of  density.  In  this  context  the  origin  of 
the  quantum  correction  to  £  is  the  same  as  the  origin  of 
the  quantum  potential. 


3.  Calculations 

The  spatial  dependence  of  s  and  the  origin  of  the 
quantum  contributions  to  transport  arise' from  gradients 
in  the  carrier  density.  These  features  are  illustrated 
through  solutions  to  the  LiouvilIc  equation  for  two 
equilibrium  solutions  using  Maxwellian  boundary  con¬ 
dition  as  discussed  in  [7].  Two  cases  are  considered.  For 
the  first  calculation  a  single  barrier  characterized  by  a 
potential 

Fo(x)=300(meV)exp[-(x/12.5A)^]  (10) 

is  placed  within  a  uniform,  ISOO  A  long,  structure  doped 
to  10*® cm”*.  The  two-dimensional  density  matrix, 
p(x,x')  as  obtained  from  the  Liouville  and  Poisson 
equations  is  displayed  in  figure  1(a).  In  equilibrium  the 
density  matrix  is  real  and  symmetric,  p{x,x')=p(x',x), 
and  the  solution  is  completely  represented  by  one-half  of 
the  matrix  on  either  side  of  the  diagonal,  x  =  x',  as 
displayed  in  figure  1(a).  The  charge  density  p(x)-p(x,x) 
is  displayed  as  a  line  plot  in  figure  1(6),  where  since  most 
of  the  structure  in  the  solution  is  contained  within  a 
range  of  250  A,  about  the  centre,  only  500  A  of  the  results 
are  displayed.  Figure  1(6)  displays  a  significant  reduction 
of  charge  within  the  barrier,  as  well  as  charge  accumula¬ 
tion  on  either  side  of  the  barrier.  While  the  reduction  of 
charge  within  the  barrier  is  a  consequence  of  the  presence 
of  the  barrier,  the  excess  charge  adjacent  to  the  barrier  is 
a  consequence  of  both  self-consistency  in  the  calculation 
and  wavefunction  (or  density  matrix)  continuity  across 
the  barrier.  The  spatial  dependence  of  the  charge  is 
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consistent  with  the  condition  of  global  charge  neutrality. 
Figure  1(6)  also  displays  two  additional  plots,  I  hc  curve 
reaching  the  lowest  value  of  density  within  the  barrier  is 
obtained  from  the  classical  Boltzmann  relation  between 
density  and  potential  energy;  the  curve  reaching  inter¬ 
mediate  values  of  density  within  the  barrier  is  obtained 
from  equation  (7).  Note  that  away  from  the  barrier  the 
density  from  equation  (7)  approaches  a  value  that  is  less 
than  the  classical  value,  a  result  that  is  a  consequence  of  a 
change  in  curvature  of  the  potential  as  the  boundaries  arc 
approached.  Neither  approximate  solution  can  be  regar¬ 
ded  as  an  adequate  representation  of  the  complete 
solution,  although  the  quantum-corrected  solution 
possesses  the  general  features  of  a  density  that  is  higher 
(than  classical)  within  the  barrier  and  lower  (than  clas¬ 
sical)  adjacent  to  the  barrier.  Figure  1(c)  displays  the 
potential  distribution.  The  lowering  of  the  potential 
adjacent  to  the  barrier  ( ~  20  meV)  is  a  consequence  of  the 
excess  charge  and  self-consistency.  Figure  1(d)  displays 
the  quantum  potential;  note  that  its  value  is  greater  than 
—  300  meV  in  the  centre  of  the  barrier.  The  energy 
density  matrix  represents  the  curvature  of  the  density 
matrix  in  the  non-local  direction.  As  seen  in  figure  1(a), 
the  curvature  is  steeper  where  there  is  excess  charge  and 
changes  sign  within  the  barrier.  The  mean  kinetic  energy 
per  particle,  £,  obtained  from  the  density  matrix  is 
displayed  in  figure  1(e),  along  with  the  Wigner  con¬ 
tribution  as  obtained  from  equation  (8).  It  is  apparent 
that  the  main  origin  of  the  structure  leading  to  the 
Wigner  contribution  is  the  quantum  potential  The 
negative  value  of  £  within  the  barrier  and  the  positive 
excess  value  of  energy  adjacent  to  the  barrier  suggest 
that  the  Wigner  contribution  is  not  a  correction,  but 
represents  a  dominant  effect,  and  that  temperature  con¬ 
cepts  (which  must  include  negative  values)  are  not  likely 
to  be  germane  within  the  context  of  equilibrium 
transport 

The  spatial  dependence  of  the  mean  kinetic  energy 
per  particle  is  also  of  significance  in  multiple  barrier 
structures.  This  is  examined  for  a  double  barrier 
structure  with 

l^(x)  =  (3(X)  mcV){cxp[  -  (x  -  75  A)/l  2. 5  A]  * 

-fexp[-(x-f75A)/12.5A]*}  (11) 

The  barriers  arc  centrally  placed  within  a  1500  A 
n*n  n*  structure  with  adjacent  10‘®cm'*n*  regions, 
and  a  centrally  placed  500  A,  10*  ’  cm  ”  *  region.  The  two- 
dimensional  density  matrix  is  displayed  in  figure  2(o). 
obtained  from  the  Liouville  and  Poisson  equations. 
There  is  excess  charge  between  the  barriers,  a  modest 
increase  in  curvature  between  the  barriers  and  a  change 
in  sign  of  the  curvature  within  the  barriers.  The  tine  plot 
of  density  is  shown  in  figure  2(6)  over  a  reduced  range  of 
6(X)  A.  The  density  as  obtained  from  equation  (7)  displays 
a  significantly  lower  charge  density  within  the  barrier  but 
order  of  magnitude  agreement  within  the  quantum  well. 
The  classical  solution  for  density  is  completely  unac¬ 
ceptable.  The  potential  distribution,  shown  in  figure  2(f), 
reaches  flat-band  beyond  400  A  on  either  side  of  the 
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Figure  1.  {a)  Density  matrix  for  a  single- barrier  structure.  The  physical  dimension  of  the  structure  is  1500  A.  requiring  that 
the  density  matrix,  which  is  calculated  over  a  square  matrix,  is  of  side  1500A/.y2.  The  centre  of  mass  and  non-local 

directions  are  indicated;  (d)  diagonal  component  of  the  density  matrix  ( - ),  from  equation  (7)  ( - ),  classical  relation 

( — _ — );  (c)  potential  energy  V{x):  (d)  quantum  potential;  (e)  energy  per  particle  from  density  matrix  ( — ),  from 
equation  (8)  (———). 


origin;  its  increase  arises  from  self-consistency  and  the 
reduction  of  charge  in  the  low-doped  region  compared 
with  the  bounding  charge  density.  The  quantum  potent¬ 
ial  displayed  in  figure  2(d)  is  positive  within  the  quantum 
well,  and  emphasizes  the  reduction  in  charge  density 
compared  with  the  classical  value;  it  is  negative  within 
the  barriers,  as  in  the  case  of  the  single-barrier  structure, 
and  positive  outside  of  the  barriers.  The  positive  value 
outside  of  the  barriers  is  a  consequence  of  wavefunction 
and  density  matrix  continuity  within  the  classically 
accessible  region.  Note  again  that  the  structure  of  the 
quantum  potential  is  apparently  the  main  origin  of  the 
structure  leading  to  the  Wigner  contribution  to  the 


energy  per  particle,  as  seen  in  figure  2(c).  As  in  the  case  of 
the  single  barrier  of  figure  1  the  calculations  suggest  that 
the  Wigner  contribution  is  not  a  correction  but  repres¬ 
ents  a  dominant  effect,  and  that  temperature  concepts  are 
not  germane  in  the  context  of  equilibrium  transport. 

4.  Conclusions 

The  calculations  of  figures  1  and  2  reveal  significant 
spatial  variations  in  energy  associated  with  density  gra¬ 
dients.  Mathematically,  these  energy  variations,  which 
are  a  consequence  of  wavefunction  continuity  as  repre- 
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Figure  2.  (a)  Density  matrix  for  a  double-barrier  structure;  (b)  diagonal  component  ot  the  density  matrix  ( — ),  from 

equation  (7)  ( - ),  classical  relation  ( - );  (c)  potential  energy  ^(x);  (cf)  quantum  potential;  (e)  energy  per  particle 

from  density  matrix  { — ).  from  equation  (8)  ( - ). 


seated  by  ‘curvature’  in  the  density  matrix,  are  suggestive 
of  quantum  heating  and/or  cooling.  Physically  these 
energy  variations  represent  the  influence  of  local 
quantum  mechanical  density  dependent  forces  on  the 
carriers.  While  their  magnitudes  indicate  that  they  must 
be  accounted  for  in  ail  quantum  mechanical  treatments 
of  transport  in  mesoscopic  structures,  an  interpretation 
in  terms  of  heating  or  cooling  in  equilibrium  is 
problematic. 
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